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The standard nonlinear perturbation theory of the gravitational instability is extended to incorporate the non- 
local bias, redshift-space distortions, and primordial non-Gaussianity. We show that local Eulerian bias is not 
generally compatible to local Lagrangian bias in nonlinear regime. The Eulerian and Lagrangian biases are 
nonlocally related order by order in the general perturbation theory. The relation between Eulerian and La- 
grangian kernels of density perturbations with biasing are derived. The effects of primordial non-Gaussianity 
and redshift-space distortions are also incorporated in our general formalism, and diagrammatic methods are 
introduced. Vertex resummations of higher-order perturbations in the presence of bias are considered. Resum- 
mations of Lagrangian bias are shown to be essential to handle biasing schemes in a general framework. 
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I. INTRODUCTION 

The large-scale structure of the universe is one of the most 
powerful probes in cosmology. The statistical nature of pri- 
mordial density fluctuations can be investigated by large-scale 
distributions of galaxies. Geometrical effects on the observed 
clustering pattern of galaxies can constrain the nature of dark 
energy, through the Alcock-Paczinski effect |l-3] or baryon 
acoustic oscillations (BAO) [4--^- Recently, it is found that 
the presence of primordial non-Gaussianity introduces scale- 
dependent bias in the halo clustering BTI-Ulll. Therefore, many 
possible models in cosmology, such as inflationary scenar- 
ios, dark energy models, modified gravity, and so on, should 
be constrained by precision measurements of the large-scale 
structure in near future. 

To compare the observations with theory, it is crucial to 
make precise predictions of observable quantities from a given 
cosmological model. The linear theory applies on very large 
scales 1 12]. However, the linear theory is not sufficiently ac- 
curate for purposes in the precision cosmology. 

Accurate predictions of statistical measures of galaxy clus- 
tering beyond the linear theory are provided by nonlinear 
theories. The method of numerical simulation is one of the 
most straightforward ways of investigating nonlinear dynam- 
ics. However, they are not free from numerical artefacts and 
systematics, such as finite-volume eff'ects, finite-resolution ef- 
fects, and so forth. Fortunately, the analytical perturbation 
theory is applicable on large scales where density fluctuations 
are small. Thus the nonlinear perturbation theory of gravita- 
tional instability attracts renewed interests in recent years. 

The nonlinear perturbation theory have been developed 
since decades ago [I3I- I20I1 . The traditional perturbation the- 
ory is formulated in Eulerian space, and such theory is re- 
ferred to as the standard perturbation theory (SPT). The per- 
turbation theory in Lagrangian space is also formulated 1I21I - 
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I29I1 . which is called the Lagrangian perturbation theory (LPT). 
The first-order LPT corresponds to the classic Zel'dovich ap- 
proximation 1.30.1 . 

In these years, the renormalized pert urbation theory (RPT) 
isil [33 l and other approaches ll33l439ll have been developed 
to improve the accuracy of perturbation theory in nonlinear 
regime, partially taking into account higher-order effects of 
the SPT. Some of those approaches are based on the reformu- 
lation of fluid equations using the propagator, the vertex, and 
a source [40], which provides a way to use standard tools of 
field theory. Nevertheless, various levels of approximations 
and ansatz should be employed in those approaches. 

The RPT and its variants mentioned above significantly im- 
prove the perturbation theory of dark matter in real space. 
However, one of the most important applications of the per- 
turbation theory is to interpret the large-scale clustering of 
galaxies or other astronomical objects, observed by redshift 
sui-veys. The observable quantity in redshift surveys is the 
distribution of objects in redshift space. Even though the RPT 
and its variants could be powerful in predicting the nonlinear 
power spectrum of dark matter in real space, one could not 
directly compare the theoretical prediction with observations. 

There are two obstacles to the comparison between the im- 
proved perturbation theories and observations. The first one 
is the redshift-space distortions: the redshift as a measure of 
the radial distance is contaminated by peculiar velocities. It is 
straightforward to take them into account in the SPT frame- 
work [141I - I44I1 . However, the SPT in redshift space breaks 
down at larger scales than in real space and the applicability 
range of scales is fairly naiTow [47], since the SPT does not 
sufficiently reproduce the nonlinear smearing effects, known 
as the Fingers-of-God (FoG) eflrect jilliell . 

Nonlinear modelings of the redshift-space distortions be- 
yond the SPT are proposed ||47| - |49| . in which the FoG ef- 
fects are phenomenologically put by hand. It is found in those 
studies that the FoG effects can be represented by putting a 
Gaussian damping factor in front of the power spectrum. It is 
shown that the Gaussian damping factor in redshift space is 
naturally derived from the LPT [i39il , where the phenomeno- 
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FIG. 1: The relation between Eulerian and Lagrangian biases. The 
Eulerian bias is expressible by the Lagrangian bias (dotted arrow) 
and vice versa (dashed arrow), only when the biases are allowed to 
be nonlocal. Note that nonlinear evolutions and formation of objects 
are nonlocal processes in nonlinear regime. 



logical Gaussian factor should be modified and additional 
mode-coupling terms should be taken into account in nonlin- 
ear redshift space. 

The second obstacle to the comparison between perturba- 
tion theories and observations is the biasing. Any astronomi- 
cal objects which can be observed are biased tracers of under- 
lying mass distributions. In the galaxy redshift surveys, the 
tracers are galaxies. The exact relationship between the distri- 
bution of mass and that of galaxies depends on the complex, 
nonlinear process of galaxy formation which is not clearly un- 
derstood. 

Analytic models of biasing have been proposed, including 
the model of local Eulerian bias ll50l - l52il . the halo model [sH- 
I57I1 . peaks model ll58l - l64ll . etc. The first model is a purely 
phenomenological parameterization of the bias, assuming a 
local relation between the mass and galaxy distributions. The 
last two models are relatively more physical than the first one, 
and are categorized as the Lagrangian bias, i.e., the locations 
of the galaxy formation are specified in initial density fields. 
The location of a peak, for example, is displaced by the dy- 
namical evolution of density fluctuations. In the framework 
of SPT, the local Eulerian bias is usually adopted |43, 44, 65- 
l67ll . However, as shown below in this paper, the halo model 
and the peaks model are not compatible with the local Eule- 
rian bias in nonlinear regime, because the gravitational evolu- 
tions are nonlocal process in general. 

In Fig.[Tl the relation between the Eulerian and Lagrangian 
biases is shown. The final mass density field is the result of 
nonlinear evolutions from the initial mass density field. The 
final number density field of objects, which are observable 
such as galaxies, is the result of nonlinear evolutions and for- 
mations of those objects. The sites of the formation in La- 
grangian space define the initial number density field. Thus 
the initial number density field also depends on the redshift of 
observation. The Lagrangian bias corresponds to the relation 
between the initial mass density field and the initial number 
density field, while the Eulerian bias corresponds to the rela- 
tion between the final mass density field and the final number 



density field. 

The initial number density field is constructed only when 
the formation process is known. For example, the halo model 
and the peaks model give prescriptions of constructing the ini- 
tial number density field from the initial mass density field; 
i.e., these models specify the Lagrangian bias. In the halo 
model with the peak-background split |53], the initial num- 
ber density is locally determined by initial mass density field. 
Thus the halo model corresponds to a local Lagrangian bias. 
In the peaks model, the sites of the formation is not locally 
determined, since the peaks are defined not only by local val- 
ues, but also by spatial derivatives of the field 1160(1 . Thus 
the peaks model corresponds to a (semi-)nonlocal Lagrangian 
bias. Such models of Lagrangian bias are more physically 
motivated than the Eulerian bias. The bias depends on the 
properties of observed objects, such as the mass and luminos- 
ity. The Lagrangian bias schemes described above actually 
depend on the mass (or the peak hight) of collapsed objects. 

The local Lagrangian bias and the local Eulerian bias are 
not compatible to each other in nonlinear regime. This fact 
is well illustrated in Fig. [1] The Eulerian bias is expressible 
by the Lagrangian bias (dashed arrow) and vice versa (dotted 
arrow). The relations involves nonlinear evolutions of mass, 
and nonlinear evolutions and formation of observed objects. 
Nonlinear evolutions and formations are nonlocal processes. 
Therefore, the Eulerian bias should be nonlocal even when 
the Lagrangian bias is local. The local biasing schemes are 
compatible only in the case that the linear theory and a local 
approximation of formation process are valid. Such condi- 
tions do not apply in generally nonlinear regime. 

In reality, the bias is definitely nonlocal. For example, the 
galaxies are largely affected by their environment [68]. The 
nonlocality should be taken into account in a precise modeling 
of the bias. The exact forms of the nonlocality in formation 
processes of observed objects have not been fully understood. 
It requires a lot of future work with analytic, numerical, and 
observational studies of complicated astrophysical processes 
to understand the exact nonlocality of bias. 

One of the main purposes of this paper is to formulate a 
nonlinear perturbation theory which can handle both the Eule- 
rian and Lagrangian nonlocal biases in general. There mainly 
two merits in this formulation. First, one can distinguish gen- 
eral properties of nonlinear clustering which do not depend 
on details of formation processes from those which largely 
depends on models of bias. Second, it is straightforward to 
predict observable quantities in any models of bias. Effects 
of redshift-space distortions and primordial non-Gaussianity 
are also included in the formalism. Resummation techniques 
in the presence of nonlocal bias are introduced as well. The 
resummations of bias are shown to be essential to handle the 
nonlocal bias in a general way. 

This paper is organized as follows. In Sec. II, both Eu- 
lerian and Lagrangian perturbation theories in real space are 
extended to include both the nonlocal bias and the primor- 
dial non-Gaussianity. Diagrammatic methods with graphical 
representations are introduced. The relation between the per- 
turbative kernels of EPT and LPT with nonlocal bias is de- 
rived as well. In Sec. Ill, the formalism of the previous sec- 
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tion is extended to include the effect of redshift-space distor- 
tions. Techniques of the vertex resummations in our formal- 
ism are introduced in Sec. IV, and finally, some models of the 
Lagrangian bias are considered, illustrating how higher-order 
bias factors are evaluated in our formahsm. 



II. NONLINEAR PERTURBATION THEORY WITH 
NONLOCAL BIAS 

A. Nonlocal bias in Eulerian space 

In the following, the density contrast of observed objects X 
at a comoving, Eulerian position x is denoted by 6-^{x). The 
observed objects X can be any astronomical objects such as 
galaxies, quasars, absorption lines, 21cm emissions, and so 
forth, which are selected and catalogued in a given redshift 
survey. The density field of the objects is not generally a lo- 
cal nor linear function of the underlying mass density contrast 
6^(x). Instead, they are nonlocally and nonlinearly related to 
each other in general. In other words, the density contrast (5x 
of objects is a functional of the mass density contrast 6^. 

The functional relation between density fields of mass and 
objects is deterministic on scales we are interested in. One 
might think that the formation process of objects is determined 
not only by the density field of mass, but also by other phys- 
ical factors such as the local radiation density and its spec- 
trum, merger histories of galaxies, etc. However, the dynam- 
ical evolutions in the structure formation are deterministic, 
and the initial density field uniquely determines all the subse- 
quent states of the universe, including the above complex fac- 
tors. Consequently, when only the growing mode solutions are 
kept, density contrasts of mass and objects are nonlocal func- 
tional of initial density contrast 8\j. we have 5m - 'FrA.^\\ 
and 5x = ;Fx[(Jl], where and Ty^ represent nonlocal func- 
tionals. 

On scales where the perturbation theory is applicable, the 
motion of dark matter is single streaming. In which case, the 
spatial distribution of dark matter uniquely inverted to give the 
initial density field 5l = I^m't^m]- Thus we have a determin- 
istic functional relation 5x[5m] - 'FyXT^\6m^. Taking into 
account the translational invariance, the Taylor expansion of 
the functional is generally given by 



5xW 



n=0 



x\ • • • Xfj bft(x — Xi , . . . , X — Xn) 



X6tn(Xl) - ■ - StniXn)- (1) 



The nonlocal bias functions, {b„}, specify the relation between 
the number density field of objects and mass density field. 
The first term with n = in Eq. ([T]i gives just a constant 
bo, which is determined by other functions b„ to ensure a 
condition {6x} - 0. The constant term bo is irrelevant as 
we consider connected moments of the density contrast. In 
Fourier space, moreover, the constant term just disappears 
when nonzero modes ^ are considered. Therefore, we 
do not retain the constant term bo in the below. 



The Taylor expansion of Eq. ([T]l is applicable only when 
the functional dependence is smooth. However, as shown in 
Sec. |IV] below, the technique of vertex resummation relaxes 
this constraint, and even singular dependences of the number 
density field on the mass density field can be handled with this 
technique. 

The local biasing ansatz is recovered when b„{x\, . . . , jc„) is 
replaced by Wr(x\) ■ ■ ■ Wr(jc„) b„, where b„ is now a constant 
for each n, and Wr{x) is a smoothing kernel with a smoothing 
radius of R. However, the Eulerian local biasing scheme is not 
a natural model when we consider nonlinear dynamics, as we 
explicitly show below in this paper 

In Fourier space, Eq. ([T) reduces to 



(2Kf6l(ki...„ - k) 



X b„{ku. . . , k„) d^iki) ■ ■ ■ 6m(k„), (2) 
where we use a notation. 



^ l- -n 



= k\ + ■ 



(3) 



throughout this paper Some variables like d^, Sx and b„ in 
Fourier space are denoted by the same symbols as those in 
real space, instead of properly using symbols like 6^, ^x, b„, 
etc. We will work in Fourier space in most of this paper. The 
convention of Fourier transform and its inverse in this paper is 
given by 



Fik) 



d'xe-'''''F{x), Fix) 



r d^k 

"J (2^ 



e*^F(jt). (4) 



Since the bias relations should not depend on the coordi- 
nates system, the bias functions b„ should be rotationally in- 
variant. For example, the first-order bias function bi{k) is ac- 
tually a function of the magnitude k = \k\ and can be denoted 
as b](k). Similarly, the second-order bias function b2{k\,k2) 
is actually a function of k\, k2, and ki2 - \ki + ^2! which 
characterize relative configuration of ki and ^2- Thus the 
function can be denoted as ^'2(^1, ^2; ^12), where the first two 
arguments are symmetric under a permutation, but the last 
one. Similarly, the third-order bias function can be denoted 
as bjiki, k2, ky, ki2, k23, k^i), and higher-order bias functions 
depend only on rotationally invariant set of variables. 

In the case of local biasing ansatz, b„{ki, . . .,k„)is replaced 
by W{kiR) ■ ■ ■ W(k„R) b„, where b„ is now a constant for each 
n, and W{kR) is a smoothing window function. When the 
smoothing radius R is much smaller than the clustering scales 
we are interested in, the smoothing window function can be 
dropped and each bias function is simply considered as a con- 
stant b„ in Fourier space. 

The mass density contrast d^nix) is also a nonlocal and non- 
linear functional of a Unear density field 6l{x). We have a 
Taylor expansion which has a similar form with Eq. (|2) in 
Fourier space: 



^(2nf6l{ki...n-k) 

izny 



xFn(ki,...,kMki)---6dkn), (5) 
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where F„'s are perturbative kernels. We consider the density 
fields at any given redshift, and time-dependences are sup- 
pressed in the above notations. For the linear density field, 
di^ik) - D(z)6o(k), where D{z) is the linear growth factor at a 
redshift z, and S{){k) is the linear density contrast at the present 
time z = 0. We adopt a normalization D(z = 0) = 1 in this 
paper. 

The evolution of the mass density field can be evaluated 
perturbatively in quasi-nonlinear regime. The SPT evaluates 
the perturbative kernel F„ in Eulerian space order by order 
|[l5,J6, 18, 20]. Forn = 1,2, for example, we have 

Fi(k) =1, (6) 

^„ , , 10 (ki h\ki-k2 4(ki-k2\^ 

(7) 

Expressions of and F4 are explicitly given in Ref. lfT6ll . 
Although above kernels F„ (n > 2) are exact only in the 
Einstein-de Sitter universe, Qm - 1, £3 a - 0, those are good 
approximations in other cosmological models with £?m 1 , 
|20]. 

Combining Eq. (|2]i and Eq. (|5]l, we have a formal expansion 
of the form 

Sx(k) - V - — ^ ■ ■ ■ —^(2n)'6'^(ki...„ - k) 

n=l \ y 

xKn(ki,...,kMki)---6dk„), (8) 

where 

Ki(k) = bi(k), (9) 
K2(kuk2) = bi(k)F2(kuk2) + b2(kuk2) (10) 
K^{kuk2,ki) = bi(k)Fi(kuk2,k3) 

+ [b2{ki,k23)F2(k2, kj) + eye] 

+ b^{kuk2,k^), (11) 

and so forth, where k - ki...„ in each expression of Kn. 

The A^-point polyspectrum P^' of the field 5x is defined by 

{6x{ki)---6x(kM)), = (2nf6l{k,...j,)P^^Hku...,kN), 

(12) 

where (■ ■ ■ )c denotes the cumulant, which corresponds to the 
connected part of the A^-point expectation value. The 2-point 

(2) 

polyspectrum is the power spectrum Px - ■ The 3- and 
4-point poly spectra are the bispectrum Bx - P]^ and the 

trispectrum Tx - Px\ respectively. Substituting the Eq. dHJ 
into Eq. ( fT2b . one can perturbatively evaluate the polyspec- 
tra f in terms of the polyspectra f of the linear density 
contrast 5l, which is similarly defined by 

{6dki)---Sdk„)}c = (27T)X(ki-n)Pf(ku...,k„). (13) 

The linear polyspectra f are proportional to the primordial 
spectra. When the initial density field is random Gaussian, 
only the primordial power spectrum is present and higher- 
order linear polyspectra all vanish. 



k I >r 

< ^ '; <^ 6„(fci,fc2,...,fc„) 




<^ f„(fel,fc2,...,fc„) 



FIG. 2: Diagrammatic rules for the Eulerian perturbation theory in 
real space. In each vertex, k = k\ + ■ ■■ + k„ should be satisfied. The 
upper and lower rules correspond to the expansions in Eq. (B and 
lO, respectively. A dashed line should be "internal": one end of a 
dashed line should be connected to a vertex with double solid line, 
and the other end should be connected to a vertex with solid lines. 



In calculating the polyspectra, diagrammatic methods are 
quite useful. Fig. |2] shows the diagrammatic rules for the Eu- 
lerian perturbation theory in real space. The first rule corre- 
sponds to each term in the expansion of Eq. (|2]l, and the sec- 
ond rule corresponds to each term in Eq. (|5). In both vertices, 
a momentum conservation ^ = ki+- ■ ■+k„ should be satisfied, 
according to the Dirac's delta function in each corresponding 
equation. The double solid line, dashed line, and single solid 
line correspond to 6x, S^, and 6l, respectively. Since the vari- 
able 5ni is expanded according to Eq. Q, dashed lines in the 
upper rule should always be connected to the vertices of the 
lower rule, i.e., the dashed lines are "internal." 

To evaluate the A^-point polyspectra of Eq. (fT2T l with ex- 
pansions of Eqs. (|2]i and (|5]), we need cumulants of the lin- 
ear density contrast which are given by Eq. (fTlT l. This 
procedure is diagrammatically equivalent to applying rules in 
Fig. [3] The open circles with cross represent the primordial 
polyspectra. When the initial density field is random Gaus- 

(ft) 

sian, higher-order polyspectra all vanish, Pj^ - for n > 3, 
and only the upper rule in Fig. |3]is relevant. 

For the evaluations of A^-point polyspectra P*^' in Eq. ( fT2b . 
we first consider A^ vertices with double solid lines. Next 
we consider possible ways of connecting those vertices with 
dashed lines and solid lines according to the rules in Figs. |2] 
and [3] The polyspectra P^' is given by the sum of terms 
which correspond to all the possible diagrams, with appro- 
priate statistical factors which are explained at the end of this 
subsection. 

Since the A^-point polyspectra are defined by the connected 
part in Eq. (fT2] i. only connected diagrams should be taken 
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k —k 

< ® > <^ Fi{k) 




FIG. 3: Diagrammatic rules for contributions from the primordial 
polyspectra. All the free ends of solid lines in Fig.|2]should be con- 
nected to each other by these rules. When the initial density field is 
random Gaussian, the lower graph with n > 3 does not exist. In the 
lower graph, ki + ■ ■ ■ + k„ = d should be satisfied. The case of n = 2 
in the lower graph is equivalent to the upper graph. 



into account. Discarding unconnected diagrams is equiv- 
alent to taking the connected part. When there exist in- 
ternal wave vectors k'^,k'2, . . which are not uniquely deter- 
mined from external wavevectors k\,...,kM, those internal 
wavevectors should be integrated with a weight of {2nY^, 
i.e., Jd^k'Jiln)^ ■ d^kyiln)^ ■ ■ ■ . The number of internal 
wavevectors to be integrated is the same as the number of 
loops in a given diagram. 

Every terms in the perturbative expansion of the polyspec- 
tra P*^' of Eq. ( fT2] i corresponds to the diagrams constructed 
by above rules. When the hierarchical orders for the linear 
polyspectra ~ ' hold, the number of loops is 

equal to the order of Pl in a given diagram. In any case, 
the perturbative order in a given diagram is apparent from the 
number and kind of crossed circles. 

When the mixed polyspectra of different types of objects, 
such as cross power spectra Pmxik), PxiXiW^ etc., are need to 
be evaluated, we just use different vertices with corresponding 
set of bias functions. 

When the bias is not present, bi - I, h2 = = ■ ■ ■ - 
0, only one dashed line can be connected to each vertex of 
the lower rule in Fig. |2] In this case one does not need to 
consider the dashed line at all, and the diagrammatic rules of 
Figs. |2] and [3] are equivalent to the ones which were previously 
introduced in Ref. il6ll . in which Gaussian initial conditions 
are assumed. Therefore, our diagrammatic rules of Figs. 
and[3]are generalization of the previous rules to the case when 
the Eulerian nonlocal bias and primordial non-Gaussianity are 
present in general. 

The diagrammatic rule for the expansion of Eq. (|8) is shown 
in Fig. m The triangle vertex corresponds to shrinking the 
vertex in terms of the diagrammatic rules of Fig. |2] In fact, for 
n - 1, 2, 3, the shrunk vertices are diagrammatically given by 
Fig. |5] These diagrammatic representations are equivalent to 
Eqs. ©-(HI]!- 

As an example. Fig. |6] shows diagrams for the power spec- 
trum up to the one-loop order in terms of shrunk vertices. It 
is a good exercise for readers to explicitly calculate the power 
spectrum from Eqs. (jSJ and (fT2l i. without resorting diagrams, 
and compare the result with the diagrammatic representation 




FIG. 4: Shrunk vertex. The triangle represents all the possible tree 
graphs constructed by the rules in Fig. |2]of EPT. The shrunk vertex 
can also be expressed by LPT diagrams. 




FIG. 5: Shrunk vertices in Eulerian perturbation theory in real space. 

of Fig.|6] The result is given by 

Px(k) ^ [Ki(k)fPL(k) 

+ ^ r ^[K2(k',k-k')fPdk')PLi\k-k'\) 
L J {zny 

+ Ki(k)Pdk) r ^K,(k,k',-k')Pdk') 

+ Kdk) [ ^K2{k',k-k')Bdk,k' -k,-k'). (14) 

One should be careful to put a correct statistical factor in 
each diagram. When there are equivalent pieces in a diagram, 
one should put a statistical factor l/«equiv^ where Wequiv is the 
number of that equivalent pieces. For example, the second 
diagram in Fig. |6] has a pair of equivalent lines connecting 




FIG. 6: Diagrams up to one-loops for the power spectrum. 



6 



both vertices. It has a two equivalent pieces, and the result- 
ing statistical factor is 1/2! - 1/2. In another way of viewing 
the statistical factor of this diagram, each vertex has a factor 
1/2! because of the prefactor 1/n! in Eq. and there is two 
ways of connecting the solid lines from the vertices, result- 
ing in the final statistical factor of(l/2!)^x2 = 1/2, which 
agrees with the previous consideration of equivalent pieces. 
In the third diagram, two solid lines in the loop connected to 
the left vertex are equivalent pieces, and the statistical factor 
is 1/2! = 1/2. In another way of viewing, left vertex has a 
factor 1 /3!, and there are three ways of choosing which solid 
lines are associated to the loop, resulting in the final statisti- 
cal factor of 1/3! X 3 = 1/2, which agrees with the previous 
consideration of equivalent pieces. The fourth diagram gives 
the same expression with the third, due to the parity symme- 
try. Similarly, the fifth and sixth diagrams have the statistical 
factor of 1/2! = 1/2 in both ways of viewing. After some 
experience, one can put a correct statistical factor in a given 
diagram. Two ways of counting as in the above examples are 
helpful for cross-checking. 



B. Nonlocal bias in Lagrangian space 

In the Lagrangian view, the dynamical evolution of cosmo- 
logical density fields is tracked by a set of trajectories of mass 
element, x{q, t), where q is the initial Lagrangian coordinates 
of the trajectory. A displacement field 'Piq, t) is defined by 



x{q, t)^q+ W{q, f). 



(15) 



and is considered as a fundamental quantity in the Lagrangian 
view of perturbations. 

Since the initial density field is sufficiently uniform, the Eu- 
lerian mass density field Pm{x) satisfies the continuity relation. 



(16) 



where is the comoving mean density of mass. On the other 
hand, the fluid elements in which observed objects reside are 
not uniformly distributed in Lagrangian space. The continuity 
relation is given by 

pAx)d^x = p\{q)cfq, (17) 

where p\{q) is the density field of the observed objects in La- 
grangian space. Note that both px and depend on the time 
of observation, because objects are identified by observers at 
a given time. The comoving mean density of the objects, px is 
common to both density fields, and the Eq. ( fTTj l is equivalent 
to the following equation: 

1 + ^x(-^) = J cPq [l + 5\{q)\ 6l[x -q- V{q)]- (18) 

When the density field is not biased, 6^{q) = everywhere. 

The expansion of the number density field by the linear den- 
sity field in Lagrangian space is formally given by 



«=i 



where b); is the n-th order nonlocal bias function in La- 
grangian space. The bias functions are essentially infinite- 
dimensional Taylor coefficients, and given by functional 
derivatives: 



,A:„) = (2 



J (^ 



)3 66dki)---S6L(k„) 

(20) 

Applying the Fourier transform to the Eq. (ITSl , and expanding 
the exponent of the displacement field, we have 

6x(k) = J d'qe-"^-^ [l + 5^(^)] - (2K)X(k) 

n\m\ J {In) 



d^k„ d'k 



dXr 
(2;r)3 

k„) 



x(2n)%(ki...„ + kl.„,-k)bl;iki, 
X ddki) ■ ■ ■ 6dkn)[k ■ nk[)] ---[k- nk'Jl (21) 

where we define = I above just for n = 0, and ¥' is the 
Fourier transform of the displacement field. In the following 
we use the displacement field both in configuration space and 
Fourier space, and we notationally distinguish between 'F and 
^. 

The displacement field 'F is similarly expanded in La- 
grangian space: 



X L„(ki 



,k„)6dki)---6M. (22) 



The evolution of the displacement field can be perturbatively 
evaluated in quasi-nonlinear regime. The LPT evaluates the 
perturbative kernel L„ order by order iIztI l28l l39l . For n = 
1, 2, we have 



L2ikuk2) 



ki ■ k2 

k\k2 



(23) 
(24) 



xb);{ku...,kn)5i^{k,)---6i^{kn), (19) 



where k - ki+k2 for L2. As in Eq. (|7]), the above kernel L2 is 
exact only for Einstein-de-Sitter universe, and the expression 
is a good approximation in other cosmological models. The 
explicit form of the third-order kernel L3 is given in iEtIi . The 
kernel L„ is not proportional to k for n > 3, in general. 

Diagrammatic rules for Eqs. (|211 and (l22t are given in 
Fig- 13 The double solid line and single solid line correspond 
to 6x and 6l, respectively. The wavy line correspond to the 
n-th order perturbation of the displacement vector 'F, and has 
an index / which corresponds to a component of the vector 
The wavy lines also carry wavevectors, and momentum con- 
servations should be satisfied at every vertices. 

Combining Eqs. (I2TI) and (l22l) . we have a formal expansion 
which should be identical to the Eq. ([8]) in EPT The first sev- 
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fcj;(fci, . . . , fc„)feij ■ ■ ■ fc,^ 



FIG. 7: Diagrammatic rules for the Lagrangian perturbation theory 
in real space. The symbols kj and L„ , are /-components of k and L„, 
respectively. The wavy lines are internal, and also carry wavevectors. 
The momentum conservation should be satisfied at each vertex. The 
diagrammatic rules in redshift space are just given by replacements 
L„ ^ LI 



eral kernels are given by 



Kiik) 
K2{kuk2) 



k ■ Li(k) + b\(k), (25) 



k ■ L2(ki,k2) + [k ■ Li(ki)][k ■ Li(k2)] 

+ b\(ki)[k ■ Li(k2)] + b\(k2)[k ■ Liiki)] 

+ b\{ki,k2), (26) 
Ki(ki,k2, ki) - k ■ Li(ki,k2, k^) 

+ {[k-Liiki)][k-L2ik2,k3)]+cyc.] 
+ [k-Li(ki)][k-Li(k2)][k-Li(k3)] 
+ {b\(ki)[k-L2(k2,k3)]+cyc.] 
+ {b\(ki)[k ■ Li(k2)][k ■ Ldk^)] + eye.} 



fe^(A;i,A:2)[A:-Li(A:3)] +cyc. 

b'i(ki,k2,k3) 



(27) 



and so forth, where k - ki...„ for K„. Those equations are 
diagrammatically represented in Fig. |8] 

When the bias is not present, bi = I, b^ = 0, and b„ - b); - 
for « > 2. In this case, the equivalence of Eqs. (l9b-(fTTTi and 
Eqs. (l25Tl-(l27Ti indicates 

Fi(Jt) = A:-Li(Jt), (28) 

F2{ku k2)^k- L2(ku k2) + [k ■ U{k,)][k ■ L,{k2)], 

(29) 

F-i{ki,k2, k^) = k ■ L^{ki, k2, k^) 

+ [[k-Li{k,m-L2{k2,k^)]+cyc.] 

+ [k ■ Li(ki)][k ■ Li(k2)]lk ■ Liik^)], (30) 



and so forth, where k - ki...„ for F„. The above equations 
give the relations of kernels in Lagrangian and Eulerian per- 
turbation theories for dark matter evolutions. Therefore, those 
equations are identities. It is easily seen that Eqs. ( l28l l. ( |29] l 
explicitly hold for Eqs. ©, (l24l l. 



The relation between Eulerian and Lagrangian bias 
functions 



In the presence of bias, the equivalence of Eqs. (l9ll-(fTTTi and 
Eqs. (|25T|-(|27]| indicates the relation between the Eulerian and 
Lagrangian bias. First several relations are given by 



bi{k)^b\{k) + l. 



(31) 



b2(ku k2) = b\(kuk2) - F2{ku k2)b\{kn) 



+ [[k ■ Ldki)] b\(k2) + (ki ^ k2)) , (32) 
b3(ki,k2,k3) = b'^{ki,k2,k3) 

- [F2(ki,k2)b\(ki2, ki) + eye] 
+ ^k ■ Li{k{)b\{ki,k2) + cyc.j , 
+ [F2{ku k2)F2(ki2, k3) + eye. 

-Fi(kuk2,k3)]b\(k) 

- [k ■ Li{k3)F2(kuk2)b\{kn) + eye] 

+ [[k ■ L2(kuk2) + k ■ Ldki)k ■ Ldk2) 

- k ■ Liikn)F2(kuk2)]b\(k3) + eye.}, 

(33) 

and so forth, where k - k\...n for bn- 

An immediate consequence of the above formulas is that 
the biasing cannot be local simultaneously both in Eulerian 
and in Lagrangian space, since the bias parameters are all con- 
stants in local bias models. For example, when the Lagrangian 
bias is local and parameters b\, b\, ... are scale-independent 
constants, the higher-order Eulerian parameters b2,b-i, . . . are 
inevitably scale-dependent according to Eqs. (|3TTi-([33]l. The 
reason for the incompatibility of local biases is that nonlin- 
ear evolutions are nonlocal process in general, as we already 
described in Introduction. 

In spite of that, local relations between the Lagrangian halo 
bias and the local Eulerian bias parameters are known in the 
halo model [57]. These relations are derived by applying the 
spherical collapse model [53, 54J, in which the density evolu- 
tions are local process, and therefore the local biases are com- 
patible. This compatibility does not hold in generally non- 
spherical collapse L75il . 



D. Comments on the bias relations in the spherical collapse 
model 



The purpose of this subsection is to find a relation between 
the bias parameters of general perturbation theory and that of 
spherical collapse model. As mentioned above, the spherical 
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collapse model is used in the halo approach to find a local 
relation between the (Lagrangian) halo bias parameters and 
Eulerian bias parameters. 

Following the same manner of Ref. |54], one can derive 
general relations between local bias parameters in Lagrangian 
space and in Eulerian space. Such relations are derived in 
AppendixlAl Up to the third order, we have 



bi^b\+ 1, 



^ ^ 1 ^ 1323 ' 



(34) 
(35) 

(36) 



where both the Eulerian bias parameters [bn] and the La- 
grangian bias parameters {b);] are local and scale-independent 
at the same time. 

Remarkably, the local bias relations in Eqs. (I34t-(l36l for 
the spherical model can be derived from Eqs. (I311-(l32t when 
the Lagrangian bias is local and Eulerian bias functions are 
averaged over directions of wavevectors. Each bias function 
b\{ki, . . . , k„) is replaced by a constant in the local La- 
grangian bias. Angular averages of Eqs. d3ll-(l3?t in this case 
are straightforwardly calculated by using the explicit forms of 
F2, Li and L2 in Eqs. (|7J, ( |23T l and (|24|. The angular aver- 
age of is given by (^3) = 682/189 l20i|96ll . As a result, 
the angular averages of Eqs. (l3ni-(l33Tl exactly reduce to the 
right-hand sides of Eqs. ([34li-(|36]l. 

Therefore, the scale- and angular-dependences of bias are 
neglected in the bias relations of the spherical collapse model, 
which are widely used in the halo approach. It is only when 
those dependences of bias are not important that the local bias 
relations of Eqs. (I34li-(l35]l are useful. In general perturbations 
without spherical symmetry, one should use the nonlocal rela- 
tions of bias in Eqs. (|3TTi-([33Tl. 



E. Comments on the stochastic bias 

In the framework of local bias models, the deterministic 
property of bias is not viable in reality. The number density of 



observed objects is not solely determined by the local density 
of mass. The stochastic biasing scheme [69.1 is a phenomeno- 
logical model to treat the biasing as a nondeterministic pro- 
cess, in the framework of local bias models. 

One of the characteristic parameters of stochasticity is the 
correlation coefficient, defined by 



(37) 



where cr„(7?) = (SiiR)}^'^, axiR) = <<5x(^))'^^ are square 
roots of the variances of density contrasts 6m(R), 6xiR) which 
are smoothed by a radius /? at a point in Eulerian space. If the 
deterministic relation 6xiR) - b{R)6m(R) exactly holds, the 
stochasticity parameter r{R) is identically unity. Deviations of 
the stochasticity parameter from unity characterize how the 
stochasticity is important. It is also common to define the 
correlation coefficient in Fourier space. 



r(k) 



Pmxik) 



^P^{k)Px{k) 



(38) 



where Pmx(k) is the cross power spectrum of mass and ob- 
jects, Pniik) and Px{k) are the power spectra of mass and ob- 
jects, respectively. 

This approach is purely phenomenological in the sense that 
the stochasticity itself does not correspond to any fundamen- 
tal physics. Instead, the stochasticity represents our ignorance 
on the formation process of observed objects. Dynamical evo- 
lutions of density field and the formation process of observed 
objects are deterministic at the fundamental level. If we would 
know all the detailed physics of the formation process, any 
stochastic character should not appear when we properly de- 
scribe precise dependences of the number density of objects 
on physical quantities, not only a local density of mass. 

For the above reason, the parameters of stochastic 
bias should be derived from the nonlocal bias functions, 
b„{ki, . . ., k„). At linear order, the correlation coefficient in 
Fourier space r(k) is identically unity, since the bias is mul- 
tiplicative, 6xik) = bi{k)6jn(k), which is a consequence of 
the translational invariance. Even in this case, the correla- 
tion coefficient in configuration space, Eq. ( |37] |. is less than 
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unity in general when the linear bias parameter bi{k) is scale- 
dependent iellToli . 

At nonlinear orders, the correlation coefficient even in 
Fourier space becomes less than unity. It is straightforward to 
calculate the Eq. dSSl l in our framework of nonlocal biasing. 
The lowest-order contribution to 1 - r{k) is given by one-loop 
diagrams. The relevant diagrams are similar to Fig. and the 
final result simply reduces to 



in real space pxix) and that in redshift space p^(s) are related 
by a continuity relation: 

p'^(s)d^s^px(x)d^x, (43) 

Therefore, the density contrast in redshift space is given by 



4(s) = [i + 5x(jc)]y-'-i, 



(44) 



1 - r(k) = 



1 



4[biik)VPUk) 

■3;,/ 



where J = d(s)/d(x) is the Jacobian of the mapping from real 
space to redshift space. One can easily calculate the Fourier 
transform of Eq. (l44l l: 



r ^[b2(k',k-k')fPdk')PL(\k-k'\). (39) 

This result is valid even when the initial density field is non- 
Gaussian. The bispectrum contributions to one-loop power 
spectra cancel out in the combination of Eq. ( l38T l. 

In a case of local Lagrangian bias, including the halo bias, 
the Lagrangian functions are constants and the second- 
order Eulerian bias function of Eq. (|32t reduces to 



6Uk) = I d^xe-"'''[l+dxix)]e- 



ik-U-(x) 



(45) 



where u- = -v,/aH and ^ is assumed. A different expres- 
sion of the above formula is seen in the Eq. (4) of Ref. liill . 
and it can be shown by partial integration that the two expres- 
sions are actually equivalent. The above relation is applicable 
even in the fully nonlinear regime. 
In Fourier space. 



4 

biikuki) ^b\ + - 



1 - 



ki ■ k2 



(40) 



u-Sk) = -ik,e(k)/k^ 



(46) 



In the large-scale limit, A; — > 0, we have b2{k', k - k') ^ b\. 



and from Eq. ( l38l l. 



1 - r(k) 



where 0{k) is the Fourier transform of the normalized velocity 
convergence, 9 = -V ■ v/aH. Expanding the peculiar velocity 
field in Eq. ( |45] ). we have 



4[bi(kWPdk) 



f p:,[Pdk')f. (41) 4(,). y r^...^f^ 



dX d^k^ dX, 



It is interesting to note that the stochasticity emerges even in 
the large-scale limit when the second-order Lagrangian bias 
parameter b!^ is nonzero. One cannot find such kind of prop- 
erties in an original approach of the stochastic biasing, since 
the correlation coefficient is just a free parameter in the lat- 
ter The halo model actually predicts the nonzero value of the 
second-order bias parameter b^ (see Sec. IV Cl below). 

To summarize this subsection, the stochastic bias is a phe- 
nomenology which is conveniently introduced in the context 
of local bias models, and is not needed in nonlocal bias mod- 
els. Stochastic properties of the local bias are derived from 
the deterministic nonlocal bias. 



III. PERTURBATION THEORY IN REDSHIFT SPACE 
WITH NONLOCAL BIAS 

A. Nonlocal bias in Eulerian space and redshift-space 
distortions 



X i2n)X(ki-n + A:'i...„, - k)^^^^^bjk[, k'J 

X e{k^) ■ ■ ■ 0(kMk\) ■ ■ ■ 6,Jk'J, (47) 

where we define bo = 1 above just for m = 0, and fi = k ■ 
z/k, fi, = kj ■ z/kj are direction cosines of wavevectors. It 
can be shown that this equation is equivalent to the Eq. (5) 
of Ref. |44], although the apparent expressions are somewhat 
different. 

In the SPT, only growing-mode solutions in each order are 
retained, and the peculiar velocity field is consistently as- 
sumed to be irrotational. Thus, the velocity field is fully char- 
acterized by the velocity divergence [20], which is expanded 
by a linear density contrast as 



m 



(2^ 



(27r)34(A:i...„ - k) 



)3 (2n)i' 
xG„iku...,k„)6dki)---6dk„), (48) 



The comoving position x in real space and s in redshift 
space are related by 1,7 1,1 



(42) 



in the plane-parallel limit of the distant-observer approxima- 
tion, where z is the unit vector along the line of sight, and v- 
is the velocity component along £. The number density field 



where / = dlnDldlna - D/HD is the linear growth rate 
which corresponds to the logarithmic derivative of the linear 
grow th factor. The perturbative kernels G„ are given by SPT 
CmililH. Forn = 1,2, we have 



(49) 



Gdk) = 1, 



G2(ki,k2)= -+{ — + —] , , +- , , h (50) 



7 \k2 k2 kik2 1 \ hk2 
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<^ G„(fci, fe2, ■ ■ • , fe„' 



FIG. 10: Shrunk vertices in Eulerian perturbation tlieory in redsliift 
space. 



k 

-Mr 



fe,> 




^ F„{ki,k2, . . . ,kn) 



FIG. 9: Diagrammatic rules for the Eulerian perturbation theory in 
redshift space. In the first rule, k = ki + ■ ■ ■ k„ + k[ + ■ ■ ■ k'^ should be 
satisfied, and k = k] + ■ ■ ■ k„ should be satisfied in the last two rules. 
The dashed lines and dotted lines are internal. 



where S„'s are derived kernels. For n - 1, 2, we have 

S i(k) ^ b,(k) + ffi\ 
S2(ki, k2) = bi(kn)F2(ki,k2) + b2{kuki) 

2Mlf^2 



+ffi^G2(kuk2) + (my 



k,k 



+ fflk 



ffitdk2) + ^bi(ki) 
k\ k2 



(52) 



(53) 



In the case of local biasing where bi and b2 are constants, 
these kernels are equivalent to the Eqs. (11) and (12) in 
Ref. |44]. Those equations are diagrammatically represented 
inFig.[ini 



B. Nonlocal bias in Lagrangian space and redshift-space 
distortions 



where the expression G2 is exact only in Einstein-de-Sitter 
universe, and only weakly depends on time in general cos- 
mology. 

Fig. |9] shows the diagrammatic rules for the Eulerian per- 
turbation theory in redshift space. The first rule corresponds 
to the expansion of Eq. (l47T i. and the second rule corresponds 
to Eq. (l48T l. The rules are used in a similar way of those in real 
space. The third iTile is common to the rule in real space, and 
coiTesponds to Eq. (|5). The momentum conservation should 
satisfied in each vertex. The meanings of double solid line, 
dashed line and single solid line are the same as in real space. 
The dotted line corresponds to the velocity convergence 9. 

Substituting the perturbative expansions of Eqs. (|5j and 
( |48] l into Eq. ( |47] i, we have a formal series of the biased field 
in redshift space: 



x5„(A:i,...,A:„)^l(A:i)---^l(A:„), (51) 



Redshift-space distortions are naturally derived in the La- 
grangian picture. The velocity of a mass element with La- 
grangian coordinates q is given by a simple derivative of 
Eq. O: 

v{q) = ax{q) = aW{q\ (54) 

where the dot represents a partial derivative by the proper time 
t. Note that we suppress the argument of the time t in our 
variables. From Eq. (|42] |. the displacement field in redshift 
space is simply given by 

z ■ ^ 

'P'^'P+ ——z. (55) 
H 

The redshift-space distoitions are exactly linear mappin gs o f 
the displacement field, even in the nonlinear regime OQIl . 
A similar equation has been applied to the analysis of the 
Zel'dovich approximation [72]. 

In our approximation that the perturbative kernels L„ are 
independent on time, we have oc D", where is the 
configuration-space counterpart of n-th order term in Eq. (l22t . 
Therefore, the time derivative of the displacement field is sim- 
ply given by 

= (56) 
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and the displacement field of each order in redshift space is 
related to the real-space displacement via 



(57) 



which is just a linear mapping of the displacement field in 
each order This linear transformation is characterized by a 
redshift-space distortion tensor s!'"^ for each n ll39ll . which is 
defined by 



(58) 



nota- 



The Eq. dSTji reduces to Wf"' = Rf^ "FjK or in a vector 
tion, •P'^*") 

As a result, each perturbative kernel of the redshift-space 
displacement is given by a linear transformation from the real- 
space kernel, 



LI = /?<"'L„ 



For n - 1 , 2, we have 



VAk) = 



k + ffikz 

^2 ' 



L\ikuk2)^- 



3 k + Iffikz 



ki ■ k^ 

k\k2 



(59) 

(60) 
(61) 



where ^ = k\+k2fo^ All the formalism of Sec HTBl in real 
space applies with substitutions of L„ — > U„ and K„ — > 5„ in 
redshift space. Corresponding to Eqs. (|25| and (I2&I . we have 

Sy{k)^k- mk) + b\{k), (62) 
S2{kuk2) = k ■ L\{kuk2) + [k ■ L\{kMk ■ V,{,k2)] 



+ b\iki)[k ■ Llik2)] + b\ik2)[k ■ Ll(ki)] 



+ b2ikuk2). 



(63) 



and so forth, where k = k[...„ for S„. Substituting Eqs. ( l60b . 
(l6Tl i. and Eqs. dSTl i. (|32] |. the above equations are equivalent 
to Eqs. ( |52] |. (1531 ). Thus, the relation of the nonlocal biases of 
Eqs. dSTT i. (|32] |. which are derived in real space, also consis- 
tently applies in redshift space. Such consistency should hold 
for any higher-order kernels. 



IV. VERTEX RESUMMATIONS 

The formalism presented so far is a natural extension of 
the SPT, simultaneously including the nonlocal bias, redshift- 
space distortions and primordial non-Gaussianity. In recent 
years, various methods beyond the SPT are developed as men- 
tioned in Introduction. In the RPT ll3Tl[32ll . the "propagator" 
plays an important role. The concept of the propagator is ex- 
tended to the "multi-point propagator" ll73[l74ll . in which the 
original propagator is identified as the one-point propagator. 

In this section, we show that the multi-point propagators 
can be obtained by resumming the external vertices in our for- 
malism. Most of the resummation methods known so far are 



only applied to the dark matter clustering in real space. Ex- 
ceptions are the Lagrangian resummation method [39i,i75i] . in 
which local Lagrangian bias and redshift-space distortions are 
included, and the time-renormalization-group method l'76'l), 
in which the halo bias is included. Therefore, the identifica- 
tion of the multi-point propagators in our formalism is an im- 
portant step toward including all the eff'ects of nonlocal bias- 
ing, redshift-space distortions and primordial non-Gaussianity 
into the resummation methods. 



A. Multi-point propagators 

We first consider how the multi-point propagator iItsIi is 
related to our formalism. We only consider one-component 
propagators with the density sector. The n-th order propagator 
r{n'(^i, . . . , k„) of the density sector is defined by an ensem- 
ble average of the functional derivative: 



, = (27Ty-'"6Uk-ki...„)Ti:\ku . . . , k„). 

(64) 

. . . ab,.^ With density 
and velocity sectors, defined in Refs. 117 3L 17411 . is re- 
lated to our definition of density propagator by - 
n!2""(Anit/-C>)'T'j'^'^...^_Mfl, ■■■Mfl„, where Anit is the linear 
growth factor at the initial time finit. The appearance of the 
Dirac's delta function on the right-hand side of Eq. (|64] | is due 
to the translational symmetry. When the initial density field is 
random Gaussian, the multi-point propagator of Eq. (|64t cor- 
responds to the coefficient of orthogonal expansion by a series 
of generalized Wiener-Hermite functionals 

Substituting Eq. (|5]l into the left-hand side of Eq. (l64l i. we 
have 



The original multi-point propagator F^"-' 



S6L(ki)---66Lik„)l (2. 



1 rd'k\ d'k' 



„)3n 2-1 m\ J {In 



)3 (2;r)3 



X (27Ty6i,(k - ki...„ - k\..JFn^,n(ku- ..kn,k\,... k'J 

x{6ak[)---6Uk'J}. (65) 

The last factor is proportional to 6^{k'-^ _^^^) for the translational 
symmetry, and the remaining factor is further decomposed 
into products of connected polyspectra, 2 n f For exam- 
ple, 

{5x626^6^) - {6x626364)^, + {5x62)c <53(^4)c 

+ <(5i53)c <52(54)c + {SxSa)^ (6263)0 

-aD p(4) + AD^Dp(2)p(2) 
- "x23Vx2i4 ^ "\2"3Vx2^i4 

+ S° S° P^2' P^2' + S° (5° P<^^ P<^^ (66) 

-I- 0)3 024^13^24 ^ "14''23^14^23 ' 

and so forth, where 6x = 6(k'^), 62 - Sik'^), 6^234 — 
{2ny'5^{k\23n), ^1234 - P(k'x, ■■ ■, ^4), etc. Thus we have 

1 rd'k[ d'k'„ 



F^^\kx,...,k„)^y - [ — 



)3 (Iny 

X FnUkx, ...,k„,k[,...,k'J {6Uk\) ■ ■ ■ 6Uk'J) . (67) 
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This equation is equivalent to Eq. (23) of Ref. ||74|] when the 
last term is decomposed into connected parts. The graphical 
representation of Eq. (l67i is similar to Fig. 2 of Ref. OTSll . 

Let us extend the above multi-point propagator of mass 
density field to include the effects of biasing and redshift- 
space distortions. We define the n-th order propagator 



, k„) of the objects X by 



= (2n)'-'"6l(k-ki...„)r'^\ku. . . , A:,,). 



Substituting Eq. ([8]) into the above equation, and following the 
same way of obtaining the Eq. ( |67] |. we have 



r^\ku...,k„)^y - f ^ 



X Kn^miku. . . k„, k\, . . . k'J {6dk\) ■ ■ ■ Sdk'J) ■ (69) 

To include the redshift-space distortions, one simply apply the 
substitution K„+„, — > S„+,„ in the above equation. The dia- 
grammatic representation of this equation is given in Fig. (TT] 

The usage of the multi-point propagator is parallel to the 
one in Refs. ll73ll74ll . The multi -point propagator corresponds 
to the summation of all the loops which are attached to each 
external vertex. Therefore, the polyspectra P^' of Eq. ( fTSl i 
are represented by using the multi-point propagators and cor- 
responding diagrams do not have any loop which is attached 
to a single external vertex of the multi-point propagator For 
example. Fig. [12] represents the power spectrum of the object 
for the Gaussian initial conditions. The resulting power spec- 
trum is given by 



d kfi 



(2nf6lik - ki...„) 



r i27Ty 



X \r'^\ku. . . , k„)\-pdk,) ■ ■ -Pdkn)- (70) 



^'x(^) 



As described in Ref. 117311 . it is important to note that each 
term in the sum is positive and the subsequent contributions 
add constructively. 



B. Lagrangian vertex resummations 

The multi-point propagators are still difficult to be exactly 
evaluated. One of the remarkable results in the RPT is a 
derivation of the propagators in the high-^ limit ll3ll [32i p73l 
p74l] . It is not trivial whether or not the same limit can be 
calculated in the presence of bias. Investigations on this line 
will be interesting for future work, and we leave them as an 
open question. Instead of evaluating the full expression of the 
multi -point propagators, we consider partial resummations of 
external vertex in the Lagrangian representation of the pertur- 
bation theory in this section. 

As obviously seen in Fig. [TT] the multi-point propagator is 
essentially a renormalized external vertex. First we consider 



partial resummation of the external vertex given in Fig. [13] In 
this figure, the gray ellipse represents all the possible graphs 
which are attached to a single external vertex with r wavy 
lines. The graphs in the gray ellipse are not necessarily con- 
nected, and can be disconnected as illustrated in the second 
line of the figure. 

According to the diagrammatic rules of Fig. [7] and appro- 
priate statistical factors, the corresponding factor of Fig. [13] 
reduces to 



I the r=o ^ \ / \ / 



X bdk\, ■ ■ ■ , k„)ki^ ■ ■ ■ kj^^^kjj 



b];(ku...,k„)ki^ 



{e-'''^)bl;{ki,...,k„)ki,---ki„, 



where 



■nk'). 



is the displacement vector at the origin. The factor 
n(k) = (e*"^) = J d^¥e-"''^n'P), 



■■h„ 

(71) 
(72) 

(73) 



is the characteristic function of the one-point distribution of 
the displacement field, where !P(!P) is the one-point proba- 
bility function of the displacement field. This characteristic 
function is a generating function of moments of !f at a single 
point in configuration space: 



i"d"n(k) 



^kj^ ■ ■ - dk 



(74) 



The characteristic function is represented bya connected mo- 
ments by the cumulant expansion theorem OTSll : 



n{k) = exp 



= exp 



^ nl 

n=l 



(75) 



In real space, the cumulants of the last line of the above 
equation are nonzero only when n is an even number for the 
rotational symmetry, and have the form 

2"n! 



{"^h ■ ■ ■ ^J2„), = [Sj^h^hh ■ ■ ■ ^j2..-,j2,. + perm.) , 

(76) 

where the factor (2n)!/2"n! is the number of all the possible 
pairings among indices . . jin, and equals to the number of 
terms in the parenthesis. When all ji,. . . ,j2n take the same 
component, e.g., z-axis, we have 

= (|!f|2"(coS 61)2")^ = ((COS0)^") 



2nH- 1 ' 



(77) 
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FIG. 1 1 : Diagrammatic representation of tiie multi-point propagator of objects in real/redshift space. 




Pxik) 



FIG. 12: Diagrammatic representation of the power spectrum with 
multi-point propagators when the initial density field is Gaussian. 
Effects of biasing and redshift-space distortions are included in the 
multi-point propagators in our formalism. 





+ . . . + 



FIG. 13: Vertex resummation with displacement legs. 



where is the polar angle of 'F, and we used the fact that the 
variables | S'l and 6 are independent and the directional cosine 
is randomly distributed in real space. 

Comparing Eqs. ( |76] | and ( ITTT i. we have 



12« 



2n + l 



(78) 



Eq. (|76] l for each n. Thus Eq. ( fTSl l reduces to 

(-1)M2„^2„ 



n{k) = exp 



n=l 



(2n)\ 



= exp 



^ (2n+ l)(2n)! 

n=l 



(79) 

When the higher-order cumulants of displacement field are 
negligible on large scales, the above equation reduces to a 
Gaussian damping factor, 



(80) 



77(^)-exp|--<im 



This equation can also be directly confirmed by contracting 



Therefore the large-scale power spectrum is smeared by non- 
linear effects, and such smearing is important in the analysis 
of BAG IHItI]. 

The above Eq. (|79T l. however, is valid only in real space 
where the clustering is statistically isotropic. In redshift space, 
the clustering is not statistically isotropic and Eq. (fTST l should 
be evaluated with the displacement field 'F^ of Eq. (l55l l. When 
the lowest-order (Zel'dovich) approximation is valid on large 
scales, and higher-order cumulants of the displacement field is 
negligible (i.e., primordial non-Gaussianity is weak), we have 
again a Gaussian damping factor, 

n(k) - exp [l + /(/ + 2),/] <|!P|2)| , (81) 

where yu = k-/k is the direction cosine of the wavevector 
with respect to the lines of sight. This damping factor rep- 
resents both effects of nonlinear smearing and FoG in red- 
shift space |39], which are present even on large scales. The 
damping factor of FoG is similar to, but somewhat different 
from that adopted in a phenomenological modeling ll47l - l49ll . 
exp(-/^/i^A:^cr^/2), where cr^ is equal to (|!Pp)/3 at the linear 
order 

Next we consider partial resummation of the external vertex 
given in Fig. [14] As in Fig. [13] the gray ellipse represents 
all the possible graphs which are attached to a single external 
vertex with r solid lines. The graphs in the ellipse are not 
necessarily connected. The corresponding factor of Fig. [14] 
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O n(k) c!;(fci, . . . , k„) A:,, ■ ■ ■ fc,,„ 
FIG. 15; Partial resummation of the vertex. 



FIG. 14: Vertex resummation with bias legs. 



reduces to 



^r\J (2;r)3 



(Pk' , 



, k„, k I, . . . , kj) 



x{6dk\)---6dk',))ki^---ki,_ 



J (2n) 



5"6Uk') 



^c];{ku 



f \d6M---6Sdkn) 



(82) 



where we define the renormalized nonlocal bias function in 
Lagrangian space: 



.,k„)^(2n)'" J 



6"6Uk') 



cPk' 



(83) 

This expression is contrasted with Eq. (|20] |. Instead of eval- 
uating the functional derivatives at 6^ - 0, taking statistical 
averages of them gives the renormalized bias functions. 

While the original function b}; can be determined solely by 
a functional relation between biased field and linear density 
field in Lagrangian space, the renormalized function c); de- 
pends on the statistical properties of the linear (initial) density 
field. In configuration space, Eq. ( [83T l is equivalent to 



c];(q-qu...,q-q„)^l^ 



S"6],{q) 



dddqi) - ■ -ddLiqn) 



(84) 



where the translational invariance is taken into account. We 
use the same symbols for variables both in Fourier space and 
configuration space as long as the notation is obvious. When 
a model of the Lagrangian bias is provided in configuration 
space, the renormalized bias functions are evaluated by the 
above equation. 

Taylor expansions are possible only when the number den- 
sity field 6^ is a smooth functional of the linear density field 
^L- However, the renormalized bias function cj; can be evalu- 
ated even when the biased field <5^ is not a smooth functional 
of 5l and does not have a Taylor expansion, since statisti- 
cal average in Eq. d83T l is possible even when the functional 
derivative in the bracket is a singular functional. We will see 




FIG. 16: An example of graphs that are not resummed in Eq. (I85t or 

Fig.ini 



some examples below in which the bias cannot be expanded 
by a Taylor series while the renormalized bias functions are 
still well-defined. 

Putting the diagrams of Fig.[T3]and Fig.[T4]together, we de- 
fine the partial resummation of the external vertex in Fig. [15] 
The gray box in this figure represents the partial resummation 
with all the sub-graphs which are attached to an external ver- 
tex with only wavy lines and with only solid lines. The partial 
resummation of Fig. [T5]results in the factor 



nik)c];iki,...,k„)ki,---ki„, 



(85) 



Connected graphs with both wavy and solid lines attached to 
an external vertex are missed in this resummation. For ex- 
ample, graphs like Fig. [16] are not included in the factor of 
Eq. ( [85] ). To obtain the full multi-point propagators, all kinds 
of graphs like the one in Fig. [T6]should be added. For exam- 
ple, the first-order propagator is diagrammatically given 
by Fig.[T7]up to one-loop contributions. 



=Q- = =o + =0-^ + + =^?^ 



FIG. 17: First-order propagator with partially resumed vertex up to 
one-loop contributions. 
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V. SOME MODELS OF THE LAGRANGIAN BIAS 
A. Local Lagrangian bias 

In the case of the local Lagrangian bias, the linear density 
field in Lagrangian space, 6^{q) is given by a single function 
of the smoothed mass density field SRiq) at the same position, 
where is a smoothing radius. In configuration space, we 
have 



Fx{Sr)., 



(86) 



where Fx is generally a nonlinear, univariate function, and the 
smoothed mass density field in Lagrangian space is given by 

5R{q)^ jd'q'WR(\q-q'\)6dq'), (87) 

where the window function Wr is spherically symmetric. Ap- 
plying the Taylor expansion and the Fourier transform to 
Eq. ( |86] |. the Lagrangian bias functions b]; in Eq. ( l20l i reduces 
to 

b];(ki, . . . , A:„) = F^^\0)W{kiR) ■ ■ ■ W{k„R), (88) 

where F^^^ - d"Fx/d6R" is the «-th derivative of the function 
Fx, and W{kR) is a (3-dimensional) Fourier transform of the 
smoothing window function Wr{x). The renormalized bias 
function c|; of Eq. ( l83l l reduces to 



c\;(ki , . . . , A:„) = {Ff(6R)) W(kiR) ■ ■ ■ W(k„R). 



(89) 



By denoting !P((5r) as the probability distribution function of 
Sr, the first factor of the above equation is given by 

d6Rp(5R)Ff{6R) 

oo 

Xoo 
d6Rp'^"\6R)Fxi6R), (90) 
oo 

where !P(«> = d"P/d6R". 

On scales which are larger than the smoothing scale, |^,| < 
\IR{i - 1 , . . . , n), one can ignore the window function WikR) 
in Eqs. ( [88] ) and ( [89] l. In this case the Lagrangian bias func- 
tions are constants: 



,k„) 



F^iO), 



(91) 
(92) 



One notices that Eq. ( [88] ) is well-defined only when the bias 
function Fx is a smooth function, because the factor F^'(O) 
corresponds to a coefficient of the Taylor series. However, 
Eq. ( [89] ) is applicable even when the bias function is not a 
smooth function. As an illustration, we consider a threshold 
bias given by 



Fx(6r) = C0(6r - (5t) ■ 
where a constant 6t is a threshold value. 



1, 



C = {0(6r - 5,)}- 



f 



dSRnSn) 



(93) 



(94) 



and is the step function. The Taylor expansion of this func- 
tion is not well-defined since all the derivatives at the origin 
are zero, F^^\0) — (n > 1), and we have b\ - Q from 
Eq. ( [88] l. However, the expectation value in Eq. ([89l l is not 
zero in this case, and using Eq. ( [90] l we have 

{Ff(6R))^{-\rC j d6Rp'^"\6R\ (95) 

for n > 1 . Therefore, the concept of renormalized bias func- 
tions c\i extends the applicability to the case when the simple 
Taylor expansion of the bias in Eq. ( [T9] l does not work. 

In the limit 6i — » oo, the threshold bias can be considered 
as an approximation to the peaks bias in the high-peak limit. 
Lagrangian statistics in this limit have been widely studied in 
1980's l.58i,,80.-821. The methods developed in those studies 
are essentially equivalent to applying the renormalized bias in 
this work to the case of local Lagrangian bias. The use of 
the partially resummed vertex of Eq. ([85) in calculating loop- 
corrections to the power spectrum is equivalent to applying a 
recent formalism developed in Ref. [751 in a case of local La- 
grangian biasing and Gaussian initial conditions. The present 
formalism is applicable even in cases of nonlocal Lagrangian 
biasing with primordial non-Gaussianity. 



B. Multivariate Lagrangian bias 

In a local bias model, the number density field 6^ is a func- 
tion of the single variable 6r of Eq. ( [87] l. We next consider the 
case in which the number density field is a multivariate func- 
tion of variables which are convolutions of the linear density 
field. The multiple variables {a - 1,2, . . .) in Lagrangian 
space are given by 



XM) = d\'UM - q')5dq'). 



(96) 



where f/„'s are the convolution kernels. For example, if one of 
the variables Xa is the linear gravitational potential 0l in La- 
grangian space, the corresponding convolution kernel is given 
hy U - -Ga^pm/\q - q'\- If one of the variables is a deriva- 
tive of a smoothed density field, say d6R(q)/dqi, the kernel is 
U = dWR(\q - q'\)/dqi, and so forth. 

The peaks bias is described by a function of the 
smoothed density field and its derivatives up to second or- 
der, (6R,di6R,didj6R) [60]. A multivariate bias model with 
two variables, Sr and is recently considered in the context 
of scale-dependent halo bias with primordial non-Gaussianity 
1^ - 85.1 . A multivariate Eulerian bias model is also proposed 

In general, we consider the biased field d^(q) is a local 
function of a finite number of the variables Xaiq) at the same 
Lagrangian position. In configuration space, we have 

Sx^Fx(xi,X2,---)- (97) 
In this case, the Lagrangian bias function of Eq. ([T9] l reduces 
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to 



dXat ■ ■ ■ dXa„ 



Ua,(kl)---Ua„(k„), 

(98) 

where Udk) is the Fourier transform of Uaiq)- The renormal- 
ized bias function of Eq. ( l83T l reduces to 



c];(ku...,K)^ ^ / 

iy< ry.. > 



^ ^ ) Ua,(ki) ■ ■ ■ Ua„(k„). 

OXa, ■ ■ ■ OXa„ I 

(99) 

In the peaks model, for example, the variables Xa con- 
tain spatial derivatives of a smoothed density field, dSR/dq/, 
d^dR/dqidqj, and corresponding kernel windows are Uaik) - 
ikiW{kR), -kikjW(kR), respectively. When one of the vari- 
ables Xa is given by the linear gravitational potential 0l, 
the kernel window Ua in Fourier space is given by Uaik) - 



C. Bias functions from universal mass function 

One of the most popular models of biasing in nonlinear 
structure formation is provided by the halo approach [SS. -iSTII . 
which is based on the extended Press-Schechter theory 11871 - 
IgTll . The peak-background split is applied in this approach, 
and the bias of halos are considered as a local Lagrangian bias. 
The bias functions can be calculated by Eq. (l88l l or Eq. (|89] l 
in this case. In this subsection, we derive explicit expressions 
of the bias functions of halos. The derivation is similar to that 
of Ref. iItsIi . in which the Gaussian initial conditions are as- 
sumed, and unfortunately the powers of growth factor in final 
expressions are incorrect. Below we correct the expressions 
of the last reference and give a derivation which applies even 
when the initial density field is non-Gaussian in general. 

The mass of halo is related to the Lagrangian radius R 
of a spherical cell by M = 47TpR^/3, or R = [M/(1.162 x 
10'2/i-iMoi3m)]'^^/2"'Mpc, where Mq = 1.989 x 10^" kg is 
the solar mass, and is the density parameter at the present 
time. Henceforth, cr^{M) denotes the variance of density fluc- 
tuations smoothed on a mass scale M which is linearly extrap- 
olated to the present time. 

According to the Press-Schechter theory and its extensions, 
the comoving number density of halos with a mass range dM 
around M, identified at redshift z, is given by 



n(M,z)dM - ■^/MF(v)iiln V, 



(100) 



where v - Sc{z)/o-(M) is the typical amplitude of fluctuations 
that produce those halos, ddz) - zfc/£>(z), and is the criti- 
cal overdensity for spherical collapse at the redshift z. In the 
Einstein-de Sitter model, the critical overdensity is indepen- 
dent on redshift, = 3(371/2)^/3/5 - 1.686, and only weakly 
depends on cosmological parameters and redshift in general 
cosmology. Since the condition of collapse is always ex- 
pressed by the linearly extrapolated overdensity at the present 



time, the growth factor is absorbed into the critical overden- 
sity 6ciz)- The multiplicity function /mf(^) is normalized by 



r 

Jo 



dv 

/mf(v) — = 1, 



(101) 



to ensure all the mass in the universe is contained in halos in 
the hmit D{z) oo. 

In the original Press-Schechter (PS) theory, the multiphcity 
function /mf(v) is given by 



/ps(v) 



(102) 



The original PS mass function is improved by Sheth and Tor- 
men (ST) llssll to give a better fit in numerical simulations of 
CDM-type cosmologies with Gaussian initial conditions. The 
corresponding multiplicity function is given by 



/st(v) = A{p) - 

71 



1 + 



1 



^/q^ 



(103) 



where p = 0.3, q - 0.707 are numerically fitted parameters, 
and A{p) = [1 H- 7r"'/^2"T(l/2 - /?)]"' is the normalization 
factor The ST mass function is applicable only for Gaussian 
initial conditions. When the non-Gaussianity is present in the 
initial density field, the multiplicity function should have the 
con-ection factor ll92ll9l . 

In the extended PS theory, the number density of halos of 
mass M, identified at redshift z, in a region of Lagrangian 
radius TJq in which the linear overdensity extrapolated to the 
present time is do, is given by iIstIi 



n(M,z\6o,Ro)dM = ^fMF(y')dln^ 



(104) 



where 



V = ^'^^^ (TO = cr(Mo), Mo = 4KpRl/3. (105) 



1/2' 



The halo of mass M is collapsed at z, while Mq is assumed 
uncollapsed at z = 0, and thus we always have ddz) > Sq. The 
conditional number density of Eq. (I104l i represents the biasing 
for the Lagrangian number density of halos. The smoothed 
density contrast do of mass modulates the number of halos. 
The density contrast of halos in Lagrangian space is given by 



n{M,z\6o,Ro) 



- 1. 



" n(M,z) 
Since c/ In v'/t/ In V = cr'^(M)/[o-^{M) - o"^], we have 



ctHM) /mf(vO 

cr2(M) - 0-2 /mf(v) 



- 1. 



(106) 



(107) 



This relation gives the function Fxi^u) of Eq. ( 186b . where the 
smoothing radius in Eq. (|86l l correspond to R ^ Rq here. We 
assume the redshift z of halo identification is the same as the 
redshift of halo observation. The smoothed linear density field 
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in Eq. ( 1861 ) corresponds to 6r D{z)6o here, because 60 is the 
value extrapolated to the present time. 

To evaluate the bias functions, the derivatives F^^ in 
Eqs. dSSb . ( [89] l need to be derived. We consider a limit of the 
peak-background split for simplicity, and assume cP-(M) » 
ctq (However, see Ref. ll94ll for limitations of this commonly 
used method.). In this limit, we have 



1 



D"(z) \d6o 



-1 



DizMM)! /mf(v) 



(108) 



in) 

for n > 1, where /^p is the n-th derivative of the multiplicity 
function /mf- The substitution 5q = is equivalent to v' = v 
in the present limit. In the same limit, taking the statistical av- 
erage over the distribution of 6q also equivalent to substituting 
v' = V in Eq. (1108b . because the distribution function of 5q is 
highly peaked at 5o = and its variance is much smaller than 
5l(z)- Therefore, we have 



4"'(0) 



-1 



D{z)cr{M)l /mf(v) 



(109) 



in the limit of peak-background split, cr^(M) » cry. The La- 
grangian bias functions of Eqs. ( [88] l, ( [89] l reduces to 



b];(ku...,k„)^c];(ku...,k„) 



-1 



DizMM)! /mf(v)' 

(110) 



where the window function to define the background field Sq is 
dropped, assuming the large-scale limit |^,| <K l//?o. The last 
approximation is consistent with that of the peak-background 
split. The above Eq. ( IllOb is applicable even in non-Gaussian 
initial conditions, as long as effects of non-Gaussianity are 
taken into account in the multiplicity function. 

In the case of ST mass function of Eq. ( 11031 ) with Gaussian 
initial conditions, the above bias functions are given by 



2p 



1 + (qv^)P 

2p(2qv^ + 2p-l) 



3qv^ + 



1 + (qv^)P 



(111) 
(112) 



and so forth. Essentially the same expressions are derived 
in Ref. JtsI] . Unfortunately, the factor 1/D"(z) is incorrectly 
missing in Eq. (55) of Ref. i75[l . Accordingly, the factor ddz) 
which appear in Eqs. (57)-(59) of that paper should be re- 
placed by Ac - D(z)6c(z). 



D. Some properties of peaks bias 

In the peaks formalism, a density peak in Lagrangian space 
is considered as a location of structure formation ll60l . In 
the smoothed density field SR(q), the number density of peaks 
above a height v is given by 

«pk = OidR/crR - v)6liV6R)\det{VV6R)\0iA3) (113) 



where ctr - (i5^)'^^, and Aj, is the smallest eigenvalue of the 
matrix [-VV(5r],7 = -djdjSR. Thus the number density of 
peaks is a multivariate function of a scalar Sr, a vector V6r, 
and a tensor VV5r at each position. 

The number density of Eq. (II 13b is a singular function. 
Therefore the Taylor expansion cannot be applied. The un- 
renormalized bias functions b„ are not well-defined in this 
case, and it is crucial to consider the renormalized bias func- 
tions of Eq. ( |83] ). The peaks bias is one of the multivari- 
ate Lagrangian bias as described in Sec. IV B| and the renor- 
malized bias function is given by Eq. ( |99l ) where (xa) - 
(6r,V6r,W6r) is a 10-dimensional vector Since W6r is 
a symmetric tensor, only six components of djdjdR are inde- 
pendent. The corresponding kernels in Eq. (|9Ft are given by 
(Ua) = [W(kR), ikiW(kR), -kikjW(kR)] where / < 

The calculation of the coefficient {d"Fx/dxai ■ ■ ■ dxa,,) for 
general n in the peaks model is quite involved. In Ref. ll62l 
|64.1 , the correlation function and the power spectrum of peaks 
are calculated up to second order in the case of Gaussian initial 
condition. In this paper, we do not derive explicit forms of the 
coefficients. Instead, we consider formal properties of the bias 
functions derived from the rotational symmetry below. 

For n - 1 , we have 



c\{k) = W(kR) 



ddRl^'^ '\d6Rj 



dSRjji 



(114) 



where Sr^ = OiSr and Srjj = 8,8 jSr. In the above equa- 
tion, the function fx is symmetrized with respect to the off- 
diagonal derivatives 8idj6R, and partial derivatives are taken as 
if didjdR and djdidR were independent when / j. From the 
rotational symmetry, the second term in the square parenthesis 
in Eq. ( II 14b identically vanishes. The last term is proportional 
to k^, since {dFx/86Rjj) oc Sij. Thus, the scale dependence of 
the first-order bias function should have a form. 



c\ik)^W(kR)(Ai+ Bik^), 



where 



dFx 
ddR 



Bi 




(115) 



(116) 



This form is exact for peaks bias models ||64ll70ll . and higher- 
order powers of k" with n > 3 do not appear. 

It is easily understood that the first bias function of peaks 
should have the form of Eq. dl 15b . The peaks are defined 
by up to second derivatives of the smoothed field, and thus the 
first bias functions of Eqs. ( |98] l, ( |99] l involve only polynomials 
of wavevector k up to second order Since the bias function is 
rotationally invariant, only the form of Eq. ( 11 15b is allowed. 

For n - 2, the same considerations show that the bias func- 
tion should have a form, 

c^iku k2) = WikiR)W(k2R) [A2 + B2 (/t^ + kj) 

+C2ki ■ k2 + D2k\kl + E2 (ki ■ k2f] , (1 17) 
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where 



A2 = 



1^ 



C2 = 



D2 



4e 



-y 



'J 



El 



ddR^iiddRjj 
30 \ dSRjidSRjj 



30 ^ 



--y 

15 ^ 

•■J 

-y 

10 ^ 



I, J 



ddRjj^ 

fFx 
ddRj/ 



(118) 
(119) 

(120) 

(121) 

(122) 



It is again easily understood that the second bias function of 
peaks should have the form of Eq. dl 17b . since the function is 
rotationally invariant and made from polynomials of wavevec- 
tors ki and ki up to second order for each. The explicit eval- 
uations of the above coefficients are tedious. In Ref. ll64ll, 
second-order biased correlation function with a Gaussian ini- 
tial condition is analytically calculated. Similar techniques 
should be also useful in our formalism, which we leave for 
future work. 



of Eulerian bias arises both from the nonlocal Lagrangian bias 
and nonlinear evolutions. 

The formalism of the present paper provides a basic 
methodology for future applications of the perturbation the- 
ory. For example, the scale-dependent bias in the presence of 
primordial non-Gaussianity has been derived by adopting ei- 
ther halo models f?-^ or the local Eulerian bias Our 
formalism allows to calculate the scale-dependent bias in any 
models of bias in a consistent manner 19511 . The power spec- 
trum with BAO in any given models of bias can be calculated 
with our formalism. The result of applying our formalism to 
a local Lagrangian bias is equivalent to the work in Ref. II75I1 . 
More precise modeling of the bias would be required in future 
analysis of the BAO in the galaxy power spectrum to constrain 
the nature of dark energy. 

Our formalism provides a way to perturbatively calculate 
the nonlinear power spectrum, bispectrum, trispectrum, and 
other polyspectra. These polyspectra are fundamental statis- 
tics and any other statistics in the large-scale structure, such 
as the correlation functions, counts-in-cells, genus statistics, 
etc. are expressible by these polyspectra, in principle. The 
formalism developed in this paper would have an essential im- 
portance in the era of precision cosmology with the large-scale 
structure of the universe. 



VI. CONCLUSIONS 

In this paper, the standard nonlinear perturbation theory of 
the gravitational instability is extended in several directions. 
One of the main extensions is the inclusion of the nonlocal 
bias, which is a general framework of biasing. Nonlocal bi- 
ases both in Eulerian and Lagrangian spaces are formulated 
and consistently included in EPT and LPT, respectively. The 
nonlinear Eulerian and Lagrangian biases are compatible only 
in the framework of nonlocal bias. The relations among per- 
turbation kernels of EPT and LPT with nonlocal biases are 
derived. 

Effects of redshift-space distortions and primordial non- 
Gaussianity are also included in our formalism. Therefore, 
our formalism provides a complete theory to predict the ob- 
servable quantities in redshift surveys, once a model of bias 
and cosmology are given. 

The concept of vertex resummations in the presence of non- 
local bias is introduced. We show that the vertex resumma- 
tion of the bias extends the applicability of the formalism to 
the case when the bias function(al) cannot be expanded into 
a Taylor series. This extension is essential for handling, e.g., 
the threshold bias and the peaks model, in which the bias in- 
volves singular functions such as the Heaviside's step func- 
tion, Dirac's delta function, etc. Calculation of perturbative 
bias functions in our formalism is exemplified by considering 
some models of local and nonlocal models of Lagrangian bias, 
such as the threshold bias model, the multivariate bias model, 
the halo model, and the peaks model. The scale dependence of 
bias functions are straightforwardly obtained in our formalism 
once a model of nonlocal bias is given. The scale dependence 
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Appendix A: Bias Parameters in tlie Spherical Collapse Model 

In this Appendix, the relations between the Eulerian and 
Lagrangian bias parameters in the spherical collapse model 
are derived, following and generalizing the argument of 
Ref. (sitj. The Einstein-de Sitter universe is assumed in the 
following equations for simplicity. The dependences of the 
results on cosmological parameters are weak. 

The time evolution of proper radius r for a spherical mass 
shell as a function of the scale factor a is given by a parametric 
formlfH 



3 1 - cos 



n 
a 
a\ 



10 
3 

10 



5i 

1/3 



(0- sin 61)2/3 



(Al) 



(A2) 



where n and a\ are initial values of r and a, respectively, and 
6i is the initial density contrast at oi. For 6i > 0, the parameter 
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6* is a positive real number. For 6i < 0, the replacement 8 iO 
should be applied to have the parameter real and positive. 

The comoving radius is given hy R = r/a, and there- 
fore an overdensity p/p of any kind in the spherical vol- 
ume is enhanced by a factor of (Ri/R)^ = (a/aiy(rj r)^ - 
{9/2)(6 - sin 0)^/(1 - cos 6*)^. Thus, the density contrasts of 
mass (5n, and of biased object 6x in the spherical volume are 
given by 1 + 5m - (Ri/R)^ in the Umit of \6i\ <K 1, and 
1 H- 5x = (Ri/R)^(l + 6^), where 6^ is the density contrast 
in Lagrangian space. Thus, we have 



9 (e-siney 



2(1 - cos 0)3' 

l+5x = (l+5m)(l+5x). 



(A3) 
(A4) 



The form of Eq. ( IA3l l is well-known fl^]. The Eq. ( IA4l i can 
also be derived from general Eqs. (IT6^ and (IT7t in the case of 
spherical perturbations. 

The linear density contrast <5l is proportional to the scale 
factor a, and from Eq. ( IA2b . we have 



The Eq. (M^ is derived in Ref. S. Although the Eq. (|A7li is 
described in Refs. lIsillSTll . they put incorrect numbers in the 
coefficients of third- and fourth-order terms. The coefficients 
fl3 and 04 in their Eq. (A4) of Ref. |54] should be replaced by 
fl3 = 2815/3969 and 04 = -590725/916389. 

The dynamical evolutions are local in the spherical collapse 
model. In this special case, both the Eulerian and Lagrangian 
biases can be simultaneously local. We have expansions 



n! 



;i=l 



(A8) 
(A9) 



where b„ and b]^ are constant bias parameters. Putting 
Eqs. (IA4l i. (|A7t , iA9i together, we have a series expansion 
of <5x in terms of (5m- Equating the resulting coefficients with 
Eq. ( IA8b gives 



3 /9\''^ 

_ (,_sin.) 



2/3 



(A5) 



For l^ml *s 1. the relation between (5m and (5l is derived as a 
power series by Taylor expansions of Eqs. (|A3) and (|A5t with 
respect to the parameter 0. The results are 
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341 



55805 



= 5i H 6, + (5i -I- 

21 567 130977 
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5l = <5m - 2Y^i 



2815 
3969 



590725 
916839^ 



(A6) 
(A7) 



bi^b\ + 1, 
b,^b\-^-b\- 



796 
1323 
7220 



^'^-T'3-^1323 



^2 + 



476320 . 



(AlO) 
(All) 

(A12) 

(A13) 



305613 

Since the relation of Eq. ( IA7l i depends only very weakly on 
cosmological model I.96.1 . the relations of Eqs. (IA10b -( IA13l l 
also do so. 
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